A class of exactly solvable models to illustrate super symmetry and test approximation 

schemes in quantum mechanics 
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We derive the analytical eigenvalues and eigenstates of a family of potentials wells with exponential 
form (FPWEF). We provide a brief summary of the supersymmetry formahsm applied to quantum 
mechanics and illustrate it by producing from the FPWEF another class of exact solutions made 
of their isospectral partners. Interestingly, a subset of the supersymmetric partners provides a 
class of exactly solvable double well potentials. We use the exact solutions of the FPWEF to test 
the robustness and accuracy of different approximation schemes. We determine (i) the ground 
state through variational method applied to an approriate set of trial functions and (ii) the whole 
spectrum using three semiclassical quantization formula: the WKB, JWKB and its supersymmetric 
extension, the SWKB quantization formula. We comment on the importance of Maslov index and 
on the range of validity of these different semiclassical approaches. 

PACS numbers: 



Analytically solvable models in quantum mechanics are 
of pedagogical interest since they allow one to illustrate 
the abstract theoretical framework by concrete examples. 
For instance, the square well potential is very convenient 
for introducing the notion of reflection and transmission 
coefficients for the scattering states with a minimum of 
calculations, yielding some astonishing results such as the 
total quantum reflection of a low incident energy parti- 
cle interacting with a potential well. In one dimension, 
besides the square potential, there are not many scatter- 
ing potentials that are discussed in standard textbooks 
of quantum mechanics. Ref. [l| gives the interesting ex- 
ample of a step potential U (x) = Uo/{l + e"^) where one 
has to deal with subtle asymptotic conditions at infinity 
to solve the scattering problem. 

Interestingly, once an analytical solution is known, the 
supersymmetry (SUSY) techniques applied to quantum 
mechanics provide a whole family of analytical solutions 
having closely related properties. 

Furthermore, analytically solvable models provide a 
testbed for comparing approximate methods to exact so- 
lutions, and can be used to model more complex situa- 
tions. 

In this article, we first investigate in one dimension 
a family of potential wells having an exponential form 
(FPWEF). In contrast with square well potentials, the 
potentials of the FPWEF are characterized by two pa- 
rameters the depth of the potential Uq and the typical 
length of variation of the potential (see Fig. [T]). We 
solve analytically the scattering problem and determine 
the bound states of the potentials belonging to the FP- 
WEF. This family provides an interesting example of the 
role of parity symmetry and the importance of bound- 
ary conditions on the existence and the number of bound 
states. 

The second part of the article is devoted to the use of 
the supersymmetry formalism applied to quantum me- 
chanics for the potential wells of the FPWEF d [1]. In 



this way, we find a new class of exactly solvable double 
well potentials. 

The third part explores the robustness and accuracy of 
approximated schemes in quantum mechanics. We com- 
pare the exact energy of the ground state for the poten- 
tials of the FPWEF with its approximate determination 
based on variational calculations. We give explicit ex- 
amples of the importance of the appropriate set of trial 
functions to get an accurate estimate of the ground state 
along with the limitations of this method. The exact 
bound spectrum of the potential wells is also compared 
with the predictions of different semiclassical quantiza- 
tion formula: the WKB, JWKB and the supersymmet- 
ric extension of the WKB formula (SWKB). We recover 
generic conclusions on the relative range of validity of 
these different approaches. 



I. THE FPWEF 

We first solve the quantum mechanical eigenstates for 
potentials Ui{x) of exponential form defined on a; > 
with a sharp wall at a; = 0, and then for even potentials 
Uii{x) of exponential form defined on the whole real axis 
and with only soft walls: Ui{x < 0) — oo and Ui{x > 
0) = -[/oe-"^, and i7ii(a;) = -[/oe-"l^l with a > and 
Ua > since we are considering potential wells. ^ The 
potential Ui{x) and Uu{x) are represented in Fig. [TJ 

The full determination of the motion of a particle of 
mass m that experiences Uioriiix) requires the knowl- 
edge of the stationary states. These states are solution 



^ The case q < has been investigated from a mathematical point 
of view in [J . 
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FIG. 1: Potentials Uu{x) and Ui{x) of depth Uo and char- 
acteristic length . Here, Uii{x) is plotted with 5 bound 
states (a = 8.48) of energies £^4 > -E3 > £2 > -Ei > -Eo (Red), 
and Ui{x) is represented for the same parameters (Black). 
This latter potential accomodates only two bound states of 
energies E3 and Ei. 



of the Schrodinger equation 
2m dx^ 



C^iorii(a;) 



(1) 



Using the dimensionless variable X = ax and the di- 
mensionless parameters a = [8mUo/{h'^a'^)]^^'^ and b = 
[8m{—E)/{h'^a^)]^^'^, we can rewrite Eq. ([!]) in the form 



d^ 

dX2 



■ 2 ~X 

a e 



(2) 



The solution of Eq. ^ can be expressed in terms of 
Bessel functions. Indeed, using the change of variable 
y = ae~'^^^, Eq. ^ takes the form of the second-order 
differential equation obeyed by the Bessel functions: 

b'Wy)=0. (3) 
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A. The bound energies and states of Ui{x) 

The energy of the bound states of the potential well 
Ui{x) are found in the negative domain of energy E < 
{b real and positive). In this case, the solution of Eq. ([1]) 
takes the general form 



iP{X > 0) = A+Jb (ae^^/^) + A+ J_ 



.-x/2 



(4) 



where Jb are the Bessel functions of the first kind. The 
spectrum is determined by the boundary conditions: 
tpiO) = 0, and ip(x 00) = 0, from which we get the 
following set of equations 



AtJbia)+A+J_b{a)^0. (5) 
[AtJbiy^O) + A+J^b{y^O)]^0. (6) 



The Bessel functions in the vicinity of zero scale as 
>/j/(2/ — >• 0) ~ y'^. As b > 0, the divergence of J-b{v) 
when y — >■ requires to set A^ = in order to fulfill 
the condition The discrete spectrum of energy {En} 
is therefore obtained from the condition ^ and involves 
the zeros of the Bessel function for a fixed value of 
the potential characteristics (depth Uq and slope a) 



(7) 



The number of bound states is thus dictated by the 
value of the index a. The ground state energy is 
Eq = —UobQ/a^. The wave function associated with the 
eigenenergy En reads V'n(a^) = ■N'nJb„{aexp{—ax/2)). 
However, if a < Oc ~ 2.405, there is no solution in b 
for Eq. ([7]) meaning that the potential Ui does not acco- 
modate any bound state. Indeed, the sharp wall of the 
potential at a; = rules out the application of the theo- 
rem according to which there is always at least a bound 
state for a one dimensional potential 0, Q . The analysis 
of the bound states of Uu{x) (see below) enables a simple 
interpretation of the non-existence of a bound state for 
Ui{x) when a < Uc- 



B. The bound energies and states of Un{x) 

The general form of the solusion of the stationary 
Schrodinger equation for the potential Uii{x) is 

> 0) = A+Jb + A+J^b , 

i^iX < 0) = A^Jb (ae^/^'j + A^ J^b (ae^^^) • (8) 

The potential Uii{x) is even and thus commutes with 
the parity operator. As a result the eigenfunctions have 
a well-defined parity. The determination of the bound 
states is made by searching for solutions such that tplx ^ 
±oo) = which implies A J — A2 = 0, and the extra 
conditions tp{0) = for the odd solutions, and ip'{0) — 
for the even solutions. For a given potential (fixed 
value of a), the corresponding discrete spectrum is given 
by the zeros of the Bessel function (for the odd 

eigenfunction) and the zeros of its first derivative 
(for the even solution) 



and 



0, 



(9) 



with feo > ^0 > ^1 > ^1 > • • ■• The eigenstates are 
^/jn{x) = A/'„J6„(aexp(— a|a;|/2)) for the eigenenergies 
En — —U[)bn/a? and ^n{.x) = A/'„ Jj,^(aexp(— Q;|a;|/2)) for 

the eigenenergies En — —Uobn/a"^ where Afn and N'n are 
normalization factors. The subset of solutions {6n} coin- 
cides with the eigenenergies of Ui{x) since they obey the 
same boundary conditions ('0(0) = and ip{x -> -l-oo) = 
0). The extra subset {&«} results from the extra sym- 
metry of the potential Uu{x) = Uii{—x) (see Fig. [Ija)). 
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FIG. 2: (color online). Reflection probability |rp for the scat- 
tering states of the potential Uii{x) as a function of the di- 
mensionless parameters a and b. The large reflection domains 
at low energy (low j3 = b/i) is a signature of quantum re- 
flection. The periodic structure is a matter-wave Fabry-Perot 
like effect. 



lution ^ yields 
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We infer the reflection probability |rp as a function of 
the dimensionless parameters a and 6 by setting A'l = 0: 
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(10) 



The result, shown in Fig. [51 exhibits quantum reflection 
for low (3 value. This is due to the fast variation of the 
de Broglie wavelength AdB(a;) — h/mv{x) when j5 tends 
to zero: 



Max 



(AX 



dB 



\ da; 



4Tr 1 



-oo. 



(11) 



Furthermore, the reflection probability displays a peri- 
odic structure as a function of a which is caused by Fabry- 
Perot cavity like resonance effect of the matter wave be- 
tween the two "walls" of the potential well. 



II. THE SUPERSYMMETRY FORMALISM 
AND ITS APPLICATION TO THE FPWEF 



The ground state is given by the first root 6o of the even 
solutions, i.e. £'o = —Uob'^/a'^. 

In contrast with the potential Ui{x) defined on half a 
space, there is always at least one bound state for the 
symmetric potential Uu{x) The threshold Uc below 
which there is no bound state for Ui{x) can now be in- 
terpreted simply through the extended potential Uii{x). 
Indeed it corresponds to the threshold for the apparition 
of the first excited state of the potential Uu{x). The 
comparison between the spectrum of Uu{x) and Ui{x) 
therefore gives a pedagogical example on which one un- 
derstands how boundary conditions influence the exis- 
tence of at least one bound state in a one-dimensional 
potential. 



C. Scattering states of Uii{x) 



The scattering states are obtained for £^ > (& = i(3 is 
purely imaginary, /3 > 0) and discussed here only for the 
potential Uu{x). An incident plane wave coming from 
— oo gives rise to a reflected and a transmitted waves. 
The asymptotic expansion for large |x| values of the so- 



Supersymmetry (SUSY) applied to one-dimensional 
problems in quantum mechanics allows to construct a 
family of exactly solvable Hamiltonians from a given solv- 
able problem [Sj]. In this section, we provide a brief re- 
minder on this method and then we illustrate it with the 
potentials of the FPWEF. 



A short reminder 



Let us consider a given solvable Hamiltonian 



H = T+V{x), where T ■ 



fi2 d^ 

— (12) 

2m dx^ ^ ' 



is the kinetic term and V{x) the potential energy term. 
The eigenvalues and eigenfunctions for the bound states 
are denoted H\il}n) = En\'4>n), with n — 0,1,2,.... Let 
us introduce the translated potential V-{x) — V{x) — 
Eq. The corresponding Hamiltonian ~ T + V-{x) 
has the same eigenfunctions, — |V'n) as H, and its 

eigenenergies, E~ — £"„ — Eq, are translated with respect 
to those of H, and are therefore positive {E^ > 0): 



H^\ij-)^E-\i^-). 



(13) 



The ground state of H- has a zero energy, i?-|V'o ) — 
so that 



V-{x) 



2m ipoix) ' 



(14) 
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where iPq{x) is the second derivative of the ground state 
wave function, Tpo{x), with respect to the variable x. The 
hamiltonian can thus be recast in the simple form 



2m 



dx^ 



<{x) 
ipoix) 



(15) 



In this form, this Hamiltonian can be somehow factor- 
ized, that is written as = , where we have 

introduced the operators 



±- 



dx ipo{x) 



(16) 



This factorization can be viewed as a generalization of 
the one developed for the analysis of the one-dimensional 
harmonic oscillator in quantum mechanics Q- Let us 
introduce the so-called superpotential 



W{x) 



(17) 



This potential is defined over the domain V of values for 
which V{x) remains finite. It has no divergence on this 
domain since the ground state wave function ^"0(2;) has 
no nodes. The relation between the superpotential W{x) 
and V- [x) is by definition 



V-{x) = W^{x) 



-.W'{x). 



This relation suggests to introduce another potential de- 
fined by 



V+{x) = W^x) + 



-.W'{x). 



(19) 



The corresponding Hamiltonian iJ+ = T + V+{x) can 
also be simply expressed in terms of the operators A^: 
H+ = A~A+. The spectrum dTS]) of iJ_ and of H+ are 
closely related. Indeed, using the expression of H± in 
terms of the operators A^ , one readily shows that 

i/_(A+|^+))=i?+(A+|V+)), (20) 
H+{A-\i^-))^E-iA-\i^-)). (21) 

As Eq = 0, we conclude that A~\iIj~) for n 7^ is an 
eigenstate of for the eigenvalues E~ . We can there- 
fore write li^^) — A^\^^), so that £'+ — E^ . Except 
for the ground state Eq , all eigenenergies of and 
coincide: i?+ = E~j^^. Starting from a given solvable 
potential with N}, bound states, one can thus construct 
by iteration a new set of exactly solvable potentials 
having respectively iVb — 1, iVb — 2, bound states. 

Supersymmetry permits one also to relate the reflec- 
tion and transmission coefficients when the two partner 
potentials, V±, have continuous spectra. Let us assume 
for sake of simplicity that the potentials V± are defined 
over the whole real axis, and that the superpotential 




FIG. 3: (color online), (a) Potential [/i (a;) with 3 bound states 
(red lines) and its first two supersymmetric partners v!^^ + 

(18) Eo and v!^^\x) + (a = 11.75). (b) Potential Uii(x) with 3 

bound states (red lines) and its first supersymmetric partner 

~(i) 



V^'' (x) + Eq that has a double well shape (a 



obeys the boundary conditions^ W{x — )■ ±00) = 0. It 
follows that V± [x — >■ ±00) = 0. We consider an incident 
plane wave e^^^ of energy E = h^k^ /2m coming from the 
direction x — > —00. Asymptotically the scattering states 
that account for the refiected and transmitted waves read 



{k,x -00) e*''^ + {k)e' 
tp^{k,x^ +00) t^{k)e 



ikx 



(22) 
(23) 



Using Eqs. and combined with Eqs. (|22j) and 
one finds r+{k) = —r^{k) and t^{k) — t-{k), which 
implies that the partner potentials have identical reflec- 
tion and transmission probabilities (|r+(fc)p = |r_(fc)p 
and |t+(fc)|2 = \Ukr). 



B. Application to Ui{x) 

In order to use the formalism of supersymmetry, we 
introduce the potential V^^\x) — Ui{x) — Eq and de- 



^ One can readily generalized the reasonning performed here for fi- 
nite but different asymptotic limits of the superpotential W{x 
±00) = W± < 00. 
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duce from Eq. (IT71) the expression for the corresponding 
superpotential W^'^^-' 



-ax/2 



J6„(ae-W2)' 



(24) 



and potential V\^' {x) = 2{W^^\x))^ + Eq - Uiix). The 

potential V_l_^'' (x) is a smooth potential that scales as ~ 
x~'^ when x tends to zero and decays as exp(— ax) for 
ax > 1 (see Fig. Ela)).^ 

Similarly, one can introduce V"i^' (x) = V"_[^^ (x) — {Ei — 

Eq) = 2(W'(i)(x))2-C/i(x) + 2£'o-£;i whose ground state 
wave function is 



(X A ipi (X 



(25) 



from which we deduce the explicit form of the superpo- 
tential: 



iy(2) = 



/2m Vo(V'oV'i - V'oV'i) 



(26) 

As previously, the potential vj^^x) = 2{W'-'^^x)f - 

(2) (2) 

V_ (x) has the same spectrum as V_ (x) except for 
the ground state. In Fig. [3ja), we represent Ui{x) for 
a = 11.75, a = 1 that accomodates three bound states 
(Nb = 3). We have also represented the supersymmetric 
partners V_1_^\x)+Eq whose two bound states correspond 

(2) 

to the first two excited states of Ui{x) and V_l (x) + Ei 
whose unique bound state corresponds to the second ex- 
cited state of Ui{x). 



C. Application to Uu{x) 

The supersymmetric potential associated with U (x) 
is directly deduced from the ground state wavefunction 
ipoix): 



^ , , {ae' 
W^^\x) = sgn(x)yC^e-"l"l/2^ 



1/2^ 



1/21 



(27) 



potential has a differentiability class C"~^^. Repeating 
the same procedure, one readily derivates the family of 
supersymmetric potential partners { } of U{x) (see 
an example on Fig. [Hb)). The supersymmetric part- 
ner vj_^'^ has a double well shape. The supersymmetry 
applied to the potential Uii(x) thus generates a family 
of exactly solvable weak double well potentials. Such a 
connection between single well and double well supersym- 
metric partners is discussed in p^ . There are not many 
examples of analytically solvable double well potentials. 
Let us mention for instance the potentials of the form 
F(x) = fc(|x| - a)2 [13. 

Note that the scattering reflection and transmission 
probabilities for the potential vj^\x) are the same as 
those of the potential Uu{x). In particular, formula ([TU]) 
gives the reflection probability for any value of the trap 
parameters encompassed in the parameter a. 



III. APPROXIMATED TECHNIQUES 

So far the results obtained are exact. In the following, 
we propose (i) to approximate the ground state energy 
of Ui{x) and Uu{x) using the variational method and to 
check and discuss the accuracy of this method and (ii) to 
test the accuracy of different semiclassical quantization 
formula for the whole spectrum. 



A. The variational method 

In order to implement the variational method, we have 
to choose a set of trial wave functions {ipcr{x)}. We will 
consider a family of trial wave functions that depend on 
a unique parameter a. The minimum of the expectation 
value of the hamiltonian H — /2m + Uioi-u{x) for these 
trial functions provides an upper bound of the ground 
state energy. The functional form of the trial ground 
state has to be chosen appropriately to get a good ap- 
proximation of the ground state energy Eq: 

Mm{{^,\H\^,)) > Eq. 



Interestingly, we can derive here the whole supersym- 
metric family as in the previous example starting from 
a potential that has a singularity in its first derivative.^ 
From Eq. ([T7]), one indeed observes that if the poten- 
tial has a differentiability class C" the supersymmetric 



^ Note that this latter potential whose analytical solution is known 
can be used for solving the 3D scattering problem for the spher- 
ically symmetric potential Ui{r) in the s-wave regime since the 
reduced radial wave function u{r) = r^{r) obeys the standard 
ID Schrodinger equation. 

* The SUSY literature contains another example of potential sin- 
gularities that can be handled: the Dirac singularities @, 0. 



1. Application to Ui{x) 

As the potential has an infinite repulsive barrier at 
X = 0, the ground state wave function vanishes at x = 0. 
First, we choose the following family of functions: 



(ficrix) 



(28) 



which is normalized to unity and obeys the same bound- 
ary condition (pa{0) = and ipa-{+oo) = as the ground 
state. This guess is inspired by the first excited wave 
function of the one-dimensional harmonic oscillator. We 
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small trap depth as illustrated in Fig. U] (green squares). 
Furthermore, the variational method allows an approxi- 
mate determination of the threshold value a^^'^^.tz ^j^g 
dimcnsionless parameter a below which there is no bound 
states. We find a™*'''''' ~ 2.5f42 as the lowest bound of a 
above which a solution of the equivalent of Eqs. pip for 
the familiy {(fcr{x)} exists. This value differs by about 
~4.5 % from the exact value. When the depth of the 
potential increases, the wave function becomes more lo- 
calized and the estimate for the ground state energy is 
much better with the trial wave functions of the family 
{(pa{x)} as it clearly appears in Fig. S) 



2. Application to Uii{x) 



FIG. 4: (color online). Relative error, A = l-Bexact — 
^5i(?7o)|/|-Bexact|, on the estimate of the ground state energy 
using the Gaussian ansatz (|32|l family {fa} (red triangles), an 
exponential ansatz {fa} (green square) as a function of the di- 
mensionless parameter a. Below the critical value Oc — 2.405, 
the one-dimensional potential Ui{x) does not accomodate any 
bound states. At the crossing of the green (square) and red 
curve (triangle), the trap is still shallow and accomodates only 
two bound states. 



now have to calculate the expectation value of H for the 
wave functions {ipa-{x)}. This quantity is a function of Uq 
and of the two dimcnsionless variables a and rj — aa/V^: 



(29) 



One finds 
Eciv) 



2m 



dtpa 



dx 



dx ~ 



3Un 



Epiv) = -Uo 



2m 4ct2 2a'^r]'^ ' 
2r] 



e'' (1 + 27/^)erfc(77) ^ 



(30) 



The minimization of the total energy is obtained for tjq 
which depends only on the dimcnsionless parameter a: 



dEi 
drj 



0, with 



Vo 



dry^ 



> 0. 



(31) 



Figure S] compares the relative error A — |-Ecxact — 
-^i('7o)|/l-£'exact| in the estimate of the ground state en- 
ergy as a function of the dimcnsionless parameter a (red 
triangles) . The poor accuracy for a small potential depth 
is due to the inappropriate functional form of the trial 
function that does not reproduce well the large exten- 
sion of the wave function in the very shallow trap limit. 
This can be confirmed by using the family of normal- 
ized wave functions of the form ifa{x) = 2xe.~^l'' j a^^"^ 
for the energy minimization. Such wave functions have a 
longer tail for large x than those of the family {Lpcix)} 
and yield a better estimate in the low a regime i.e. for a 



The minimization of the hamiltonian expectation value 
for the potential \J\\(x) is here performed in the sub-space 
of Gaussian trial functions 



(p.{x) = ( ^ 



1/4 



1 



rl/2 



(32) 



This family of trial functions is inspired by the ground 
state wave function of the one-dimensional harmonic os- 
cillator and has no node as expected for the ground state 
of a potential well. We find, for example, E{r]o) ~ 
— 0.545?7o for a = 5 which differs from the exact value 
by about 1 %. 

The variational method also allows for the determina- 
tion of first excitated state. For this purpose, one has to 
choose a family that has the same symmetry as the state 
considered and that is orthogonal to the family of states 
used for the determination of the ground state. From 
this respect, the extension of the trial functions used for 
the potential Ui{x) to the family { <Pa-ix > 0) — '^a{x) 
and '^(t{x < 0) = —ifrj{—x) } provides a new family that 
is precisely orthogonal to the family { (^^r }, has one node 
and is thus appropriate for the determination of the first 
excited state through the variational method. This cal- 
culation is exactly the one performed on the half space 
a; > for the potential Ui(x) whose accuracy is sum- 
marized on an example in Fig. U) Thereby, the search 
for an approximate of the ground state energy of the po- 
tential Ui{x) through the variational principle gives an 
estimate of the first excited state of the symmetric ex- 
tension Uii{x) of the potential Ui{x). 



B. The semiclassical quantization formula 

In order to get an approximate determination of the 
whole spectrum, one relies on semiclassical quantization 
formula. In this section, we test the accuracy of three 
semiclassical approaches. 

The most commonly used is the WKB quantization 
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FIG. 5: (color online). Comparison of WKB (blue squares), 
SWKB (green disks) and JWKB (red diamonds) predictions 
for the energy spectrum of Ui(x) with a — 32 (10 bound 
states) with exact results. The relative error A is plotted 
for the 10 bound states as a function of the dimensionless 
parameter a. 

condition which reads^ [l2l - [T5j | 

j> pdx= {n + u)h, (33) 

where v is the Maslov index, whose value is determined 
by the matching conditions of the wave function at the 
turning points. A smooth (sharp) wall gives a contri- 
bution 1/4 (1/2). For the potential Ui{x), we have a 
sharp wall at x = and a smooth one for x > so that 
i/ = 1/4 + 1/2 = 3/4. The calculation of the action § pdx 
in the case of the potential Ui{x) combined with Eq. ([55)1 
gives the following implicit equation for the approximate 
semiclassical determination of the eigenenergies: 

+ J-nyn), (34) 

with F{y) = yjl-y"^ - ycos'^(y), = exp(-ax„/2) 
and U{xn) = E. Semiclassical approaches are supposed 
to work better for large quantum numbers. We compare 
in Fig. [5] the exact energies of the potential Ui{x) with 
a choice of parameters (a = 32 and a — 1) such that 
it accomodates 10 bound states with the approximated 
values obtained from Eq. We indeed observe an 

accuracy which gets better and better up to the fifth 
level but then which get worse. 

The standard WKB estimate for the energies is sig- 
nificantly improved by taking into account higher orders 
corrections in h to the standard WKB quantization con- 
dition, this approximation scheme is referred to as the 



JWKB quantization condition [T^. The first correction, 
(5, reads 

[n + v) = ^<j> pdx + 5 (35) 

where 

The explicit form for the potential Ui{x) reads: 

Compared to the WKB results, we obtain an improved 
accuracy for the whole spectrum (see Fig. [S]). 

Combining the supersymmetry formalism with the 
WKB method, one can work out a SWKB quantization 
condition [TtI . [l8l | . This third semiclassical quantization 
formula reads 

/ [2m{E- - W^ix))] ^''^ dx = nhn (37) 

where E' = W^{xrain) = W^^(a;max)- The SWKB ap- 
proach yields even the exact bound state spectra for all 
shape invariant potentials (SIPs), that is when the pair 
of the SUSY partners V+ and V- are similar in shape 
and differ only in the parameters [l^ [2^ . The potential 
Ui{x) does not belong to the SIPs and thus provides an 
interesting example to test the accuracy of the SWKB 
spectrum prediction. By construction, the SWKB ap- 
proach requires the knowledge of the ground state wave 
function and thus gives the exact ground state energy. As 
a result it provides the best estimate for the deep energy 
states as it clearly appears in Fig. [S] Among the three 
semiclassical approximation schemes, the JWKB formula 
turns out to be the most accurate for the states near the 
continuum. These conclusions on the relative range of 
validity and accuracy of the different semiclassical quan- 
tization formula is quite generic. 

Let us consider the highest bound state for a deeper 
potential.^ For the three approximation schemes (WKB, 
JWKB and SWKB), the deeper the last bound state the 
better the estimate. The JWKB quantization formula 
gives systematically a better account of the energy of 
the last state. This result is well known in molecular 
physics. Note that the JWKB quantization condition for 
the highest vibrational levels in a molecular potential can 
be further improved using the Gribakin and Flambaum 
formula [2l[ for the scattering length [22l. f23l|. 



^ For a ID harmonic potential or a Morse potential, this approxi- 
mation gives a spectrum that coincides with the exact one. 



® In 3D quantum scattering theory, this state plays a particular 
role since it determines the magnitude and sign of the scattering 
length, a result also referred to as the Levinson's theorem |21| . 
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The WKB quantization rule for the potential Uu{x) 
yields 

+ = (38) 

with Un = exp(— aa;„/2) and U{xn) is the energy of the 
n'^ state since the potential well Uu{x) has two smooth 
walls (z/ = 1/4 + 1/4 = 1/2). As expected, the odd 
values of n coincide with the energies determined by ap- 
plying the WKB quantization condition with Ui{x) (see 
Eq. ([M])). This example provides a clear illustration of 
the importance of the Maslov index. An extra subset of 
energies is obtained that corresponds to the even values 
of n, it includes the ground state n = which has an 
energy lower than the ground state Ui{x) (see Fig. [T]). 

IV. DISCUSSION 

We have shown the application of supersymmetry to 
the potentials Ui{x) and Uii{x) that were not investigated 
so far in this field, and illustrated the importance of ex- 
act solutions to test approximate methods in quantum 
mechanics. The same approach can be used to analyze 
the family of potentials of the form In addition to the 



pedagogical value of these examples, we would like to em- 
phasize that supersymmetry complements perfectly the 
traditional teaching of quantum mechanics at the under- 
graduate physics. It answers precisely to basic and im- 
portant questions such as: can two potential wells have 
the same spectrum and different shape ? Is it possible 
to rebuild the potential shape knowing its reflection and 
transmission probabilities for all incident energies ? Does 
it exist transparent potential ? It generalizes the factor- 
ization procedure of the stationary Schrodinger equation 
introduced for the harmonic oscillator, enlarges the class 
of exactly solvable potentials and provides new approxi- 
mating schemes for quantization rules. 
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